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Abstract 

The hyperfine splitting for the 2S charmonium state is calculated and the predicted number 
is Ahf(2S') = 57 ± 8 MeV, being by derivation the lower bound of this splitting. It results 
in M{rj c {2S)) = 3630 ± 8 MeV, which is smaller by two standard deviations than found in the 
Belle experiment [lj], but close to the n c (2S) mass observed by the same group in the experiment 
e+e" -» J/ipr) c 6] where M{rj c (2S)) = 3622 ± 12 MeV was found. 

PACS numbers: 11.15.Bt, 13.25. Gv, 14.40.Gx 
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I. INTRODUCTION 



Recently the Belle Collaboration has observed a new charmonium state, the r] c (2S), in 



exclusive B — > KKsK~ir + decays l|. The measured mass of the r) c (2S), M(2 1 Sq j 
3654 ± 14 MeV, has appeared to be rather close to the ip(2S) mass, M 2 = 3686 MeV, giving 
rise to a small hyperfine splitting for the 25 state, 

A HF (2S, exp) = M 2 - M(2 1 S ) = 32 ± 14 MeV. (1.1) 

Compared to the hyperfine splitting for the IS* state, Ahf(15', exp) = 117.2 ± 1.5 MeV 
the number (jl.lj) is 3.5 times smaller, 

A HF (l^exp) V ; 

and has a large experimental error coming from the error in the r] c (2S) mass. However, in 
the ratio (jl.2j) even the upper limit is rather small and it is of interest to compare this small 
ratio with the predictions in the conventional theory of the spin-spin interaction in QCD. 

In this rapid communication we shall discuss the problem posed in a model-independent 
way and show that the experimental value (|1.2j) puts a strong restriction on the coupling 
«hf(2S'), determining the spin-spin interaction. In particular, the central value of Ahf(2>S) 
corresponds to o;hf(25') = 0.18^-0.19 which implies the very large value for the renormaliza- 
tion scale fi2 = 2M2 = 7.4 GeV. This scale appears to be drastically different from that for 
the 15 state (J/ip) where the scale in the strong coupling determining the HF interaction 
is fii = \M\ = 1.55 GeV (Mi = M(J/tp)) and the coupling constant ctjuK/^i) — 0.335 is 
rather large. We shall show that if one uses for the 2S state the same prescription as for 
J/i/j, then our theoretical prediction is Ahf (25, theory) = 57 ± 8 MeV, which by derivation 
is the lower bound of A H f(2S') and gives rise to the mass value M(r] c (2S)) = 3630 ± 8 MeV, 
being two standard deviations off the value obtained from the Belle experiment 

II. MODEL-INDEPENDENT CALCULATION OF A HF (2S) 

In one-loop approximation in the MS renormalization scheme, the hyperfine splitting, 
Ahf(tiS'), is given by the well-known expression: 

a / m 8a H F(/i™) |D m|2 A . OtUF^n) c \ fr) ^ 

A HF (nS) = — \Rn{0)\ 1 H 4hf , (2.1) 

9 m A c \ 7T J 



where the factor £hf comes from the one-loop corrections, £hf = J2P0 — § — ^ In 2 [3J, and 
for rif = 4, ^hf(^/ = 4) = 0.2857 is small, so that the ckhf correction in the brackets turns 
out to be less than 3% (aHF^n) < 0.35) and can be neglected in the ratio (|1.2jl . 

For light mesons the relativistic version of the expression (12.1)1 also exists, it was derived 
in Refs. j^J and can be useful for charmonium, since the current c-quark mass is not large, 
m c = 1.3. ±0.2 GeV Q, 

KAnS) = l^^\R n (0)\ 2 . (2-2) 



Here u n is the kinetic energy matrix element: 



u n = (\fp 2 + m 2 c ) nS , (2.3) 

which plays the role of the constituent quark mass. The HF splitting (|2.1|) or ()2.2|) strongly 
depends on the c-quark mass chosen and io\ can change by a factor of 2 for different choices 
of m c . To escape the problem of a correct choice of m c it is convenient to consider the ratio 
of hyperfme splittings with relativistic corrections: 

R 2 (0) 



Ag F (2£) = «HF(/j 2 ) f^i\ 2 
Ag F (lS') «hf(/Ui) V^2/ 



(2.4) 



Ri(0) 

Note that the difference between uji and uo 2 in charmonium which takes into account rel- 
ativistic corrections, is not large, e.g., for the spinless Salpeter equation (cui/c^) 2 = 0.94, 
however, it will be useful to keep this factor in our analysis, since due to this factor the ratio 
f!2.4|) and therefore the hyperfme splitting is 6% smaller. However, the difference between 
u)(nS) and the current (pole) mass m c coming from relativistic corrections, is not small, 
being about 200 MeV and 250 MeV for the IS and 2S states, respectively (see Section ITVj) . 

The ratio of the wave functions at the origin occurring in Eq. f!2.4|) can be extracted in a 
model-independent way from the leptonic widths: 



where the factor 7„ is 



4a 2 e 2 

(nS) = —^{Rnm 2 ^, (2.5) 

n 
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7n = l-— a s (M n ), (2.6) 

D7T 



and in 7„ the strong coupling a s is taken at the scale n n = M n {Mi = M(J/i/;), 
M 2 = M(ip(2S)). Note that just this prescription for a s (/j,) (in the QCD factor 7) pro- 
vides the minimal value of the extracted wave function at the origin. For A( 4 )(3-loop) = 280 



3 



MeV (rif = 4), which corresponds to a s (M z ) = 0.117, the values of a s (M n ) in 3-loop ap- 
proximation are the following, 



a s {M x ) = 0.247, a s (M 2 ) = 0.232, 
7i = 0.581, 7 2 = 0.606. 



(2.7) 



Note that for another choice: 71 = 7 2 which is often used, the value of A H f(25') would be 
4% larger than in our consideration. 

Then from the ratio of leptonic widths (|2.5|) it follows that 



#2(0) 



M 2 \ 2 rw(2S) 7l 



i?i(o) Vm x ; r e+e - (15)72 ' 

Taking the experimental values: r e + e -(lS) = 5.26 ± 0.37 keV and T e 
keV |2j, this ratio with a good accuracy is 

#2(0) I 2 



(25) 



Ri(0) 



(0.59 ±0.08)—. 

72 



(2.8) 
2.19 ±0.15 

(2.9) 



The values of the wave functions at the origin extracted from the leptonic widths are the 



following: |i?i(0)| 2 = 0.917 GeV 3 and \Ri 



0.518 GeV 3 and at this point it is important 



to stress that for the 2S state the extracted value of the wave function at the origin already 
takes into account the influence of the close lying DD channel on the wave function in an 
implicit way, and it is the same both for the hyperfine splitting and the leptonic width. 
Due to the nearby DD threshold the experimental value of |i?(0)| 2 as well as the leptonic 
widths typically appear to be smaller by about 20% than in the theoretical calculations 
(one-channel approximation) and the influence of the DD channel on the ip(2S) mass was 
discussed in Ref. 0]. 

Now combining the expressions (|2.9|) and (J2.4|) one obtains 

Ag F (25) 



Agp(lS) 



(0.59 ± 0.08)t/hf 



with the factor 



'wA 2 a HF (Ai2)7i 
Vhf = — I 



(2.10) 



(2.11) 



V U> 2 / "HF(Atl)72 

From the analysis of the leptonic width and hyperfine splitting for the J/ip it is known 
that the renormalization scale \i\ in aHF(lS') is different from M\ and a good description of 
the Ahf(IS') can be reached if in Eq. ()2.1j) or ()2.2j) the scale is taken to be \i\ = 5 Mi = 1.55 



GeV, which with A( 4 )(3 - loop) = 280 MeV gives the value Q!hf(|^i) = 0.335. Now with 
the same prescription /i 2 = hM 2 , the value «hf(|^2) = 0.307 is obtained, and then from 
Eqs. (JHZJ) and (l2~TT|) it follows that 

r7 HF = 0.826, for (— V = 0.94, ^ = 0.96. (2.12) 

\uj 2 J 72 

As a result from Eq. ()2.10|) the "theoretical" ratio of the hyperfme splittings is 

Ag F (2S) 



and our prediction for the Ahf(2S') is 



(0.48 ±0.07) (2.13) 



A HF (25, theory) = (57 ± 8)MeV. (2.14) 

Owing to our derivation this number can be considered as the lower bound of Ahf(2S") which 
was obtained in a model-independent way. The predicted mass of the r) c (2S), M(t) c (2S)) = 
3630 ± 8 MeV, differs by two standard deviations from the measured value in the Belle 
experiment jlj], but is rather close to the t] c (2S) mass observed in another Belle experiment 
(f| e + e~ — > J/ipr] c where the measured rj c {2S) mass, M(j] c (2S)) = 3622 ± 12 MeV is in good 
agreement with our prediction. Note that besides the dominant perturbative term in the 
HF splitting there exists also a nonperturbative (NP) contribution to Anp(nS). However, 
the NP terms in charmonium can be calculated as in Ref. \\ and turn out to be small: we 
find Ag^lS) = 3-5 MeV and Ag^S) = 1-2 MeV. It is of interest to note that the 
same model-independent estimate of the Ahf(2S') only with 71 = 72 and «hf(a*i) = «hf(a*2) 
was suggested many years ago at the time when the mass of the rj c (lS) was not correctly 
measured 



III. THE RENORMALIZATION SCALE IN a nF (p) 

To obtain the central value in Eq. (jl.lj) in the same procedure one needs to take the very 
small value ohf^) — 0.187, which corresponds to the large renormalization scale: /i 2 = 7.4 
GeV > 2M 2 . It is difficult to point out any physical explanation for such different scales of 
H2 = 2M 2 for the 2S state and and \i X = \M l = 1.55 GeV for the J/if>. 

Existing theoretical calculations of Ahf mostly provide the 2S splittings in the range 
70-90 MeV Q, Q, while only in Ref. [llj . where the modified (screened) color Coulomb 



interaction was used, the calculated HF splitting is smaller: Ahf(2S') = 38 MeV. However, 
in Ref. d| with the same modified Coulomb interaction the fine structure splittings for the 
Xc(l-P) mesons turned out to be twice as small as the experimental data. 



IV. CHARMONIUM SPECTRUM 

Here we give as an illustration the spectrum and the constituent masses u n in charmonium 
which were calculated solving the spinless Salpeter equation with a linear plus Coulomb 
potential, 



the parameters a s t = 0.42 and a = 0.18 GeV 2 were taken. For the one-loop pole mass 
we used m c = 1.42 GeV. The remarkable feature of this set of parameters is that the 
meson mass, M(nL) = M^\nL) - C SE (nL), contains a very small negative subtractive 
constant which can be strictly determined by the nonperturbative self-energy contribution: 
Cse = —4(J ■ 0.2A/iru! n L [3 an d its value Cse = —33 ± 3 MeV is small and approximately 
equal for all nS states (n < 5). 

From Table |J one can see that the energy u {nS) of the c quark, playing the role of the 
constituent mass of a given nS state, appears to be around 1.65-1.75 GeV, i.e., essentially 
larger than the current mass, m c = 1.4 GeV in the Salpeter equation and gives rise to a 
small suppression of the HF splittings. This fact can be considered as an explanation why in 
the nonrelativistic approach the constituent c-quark mass is usually taken to be larger, e.g., 
in Ref. Q m c = 1.84 GeV. Note also that with u(lS) = 1.60 GeV and neglecting the NP 
contribution of about 3-5 MeV, the value of (x-hf(1S) = 0.355 is needed in order to obtain 
the experimental number for Ahf(IS') = 117 MeV. 

V. CONCLUSIONS 

Thus we can conclude that the Belle experiment on the r) c {2S) mass if it is confirmed, 
creates a number of interesting theoretical problems, in particular, about the correct choice 
of the renormalization scale of a s (/x) for the excited states. In the conventional approach 




(4.1) 



In the static potential 



v o[r) = ~- h err 

3 r 



(4.2) 
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TABLE I: Spin-averaged masses M(nL) for the spinless Salpeter equation with the static potential 
V (r) (m c = 1.42 GeV, a = 0.18 GeV, a st = 0.42). 



O) rp 





Ll?*, r ( in CtpV ) 


M(nL) a ) in MeV 
theory 


M(nL) in MeV 
exp er iment 


IS 


1.64 


3069 


3067 6±0 05 


2S 


1.69 


3662 


3678. 0±2. 6 


3S 


1.75 


4083 


4040±10 


4S 


1.80 


4433 


4415±6 


IP 


1.63 


3526 


3525.3±0.2 


2P 


1.70 


3967 




ID 


1.66 


3825 


(l 3 L»i state) 








3769.9±2.5 


2D 


1.72 


4200 


4159± 20? 



re overall self-energy constant Cse = —33 MeV is used. 

aHF^) = 0.30 and A HF (25) = (57±8) MeV and our prediction for the mass is M(r) c (2S)) = 
3630±88 MeV which appears to be in agreement with the result of the other Belle experiment 
e + e~ — > J/ip,r] c . To obtain decisive conclusions it would be important to measure the 
i] c (2S) mass with better accuracy. The confirmation of the value of A-ap(2S) in the range 
close to 32 MeV would require a drastic reconsideration of our understanding of the scale 
in the HF interaction and possibly also of the QCD factor 7„ present in the leptonic width, 
since they are interconnected. 
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